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^ ABSTRACT 

Cn Context. The southern long-period Cepheid RS Pup is surrounded by a large circumstellar dusty nebula reflecting the light from the 

^^ central star. Due to the changing luminosity of the central source, light echoes propagate into the nebula. This remarkable phenomenon 

was the subject of Paper I. The origin and physical properties of the nebula are however uncertain: it may have been created through 
■j' mass loss from the star itself, or it could be the remnant of a pre-existing interstellar cloud. 

Ph Aims. Our goal is to determine the three-dimensional structure of the light- scattering nebula, and estimate its mass. This will bring 

^0 us new clues on the origin of the nebula. Knowing the geometrical shape of the nebula will also allow us to retrieve the distance of 

RS Pup in an unambiguous manner using a model of its light echoes (in a forthcoming work). 

Methods. The scattering angle of the Cepheid light in the circumstellar nebula can be recovered from its degree of linear polarization. 
We thus observed the nebula surrounding RS Pup using the polarimetric imaging mode of the VLT/FORS instrument, and obtained a 
Q map of the degree and position angle of linear polarization. 

^ Results. From our FORS observations, we derive a three-dimensional map of the distribution of the dust around RS Pup, whose 

^ overall geometry is an irregular and geometrically thin layer. The nebula does not present a well-defined central symmetry. Using a 

^ simple scattering model, we derive a total dust mass of Mdust = 2.9 + 0.9 M© for the light- scattering dust within 1.8' of the Cepheid. 

>— ' This translates into a total mass of Mgas+dust = 290 + 120 M©, assuming a dust-to-gas ratio of Mdust/^gas = 1.0 + 0.3%. 

Conclusions. The high mass of the dusty nebula excludes that it was created by mass-loss from the star. However, the thinness of the 
^~^ dust distribution is an indication that the Cepheid participated to the shaping of the nebula, e.g. through its radiation pressure or stellar 

^ wind. RS Pup therefore appears as a regular long-period Cepheid located in an exceptionally dense interstellar environment. 

OO 

■^- Key words. Stars: individual: RS Pup, Stars: circumstellar matter. Techniques: polarimetric. Stars: variables: Cepheids, ISM: dust, 

V^ extinction. Scattering 
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^^ 1. Introduction ment. These data give us accurate maps of the polarization de- 

(N( gree and angle over the nebula (Sect.l2]i. The polarimetric prop- 

^H RS Pup (HD 68860, SAO 198944) is the only Cepheid known to grties of the dust are related to the scattering angle of the light, 

'' be mtimately associated with a large circumstellar nebula. From ^nd thus the position of the scattering dust relative to the plane of 

. ;:; NTT imaging of the light echoes present in the nebula, Kervella the sky. We complete these data by the NTT/EMMI observations 

^ et al. (EQQ8, hereafter Paper I) derived its distance with a preci- presented in Paper I (Sect. ^. We then estimate the thickness 

jH sion of ± 1 .4%. The method employed in that article relies on the of the light-scattering dust and its three-dimensional distribution 

03 assumption that the observed nebular features are located close (Sectjfli. We discuss the total mass and origin of the nebula in 

to the plane of the sky. As argued by Bond & Sparks (2009), Sect.pfin AppendixlAl we examine the possibility that the neb- 

this assumption, if not correct, could lead to a significant bias ui^r features evolved over the last five decades, 
on the determined distance. In order to overcome this limita- 
tion, the goal of the present work is to determine the three- 
dimensional geometry of RS Pup's nebula. For this purpose, we 

obtained imaging polarimetry data using the VLT/FORS2 instru- ^ VLT/F0RS2 polarimetric observations 



Send offprint requests to: Pierre Kervella 2.1. Instrumental setup and observation log 

* Based on observations made with ESO telescopes at La Silla 

Paranal Observatory, under ESO program 383. D-0623 (A). We observed RSPup on the nights of 4 April, 26 April and 

Correspondence to: pierre.kervella@obspm.fr 26 May 2009 using the FOcal Reducer and low dispersion 
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Table 1. Log of the FORS observations of RS Pup. 






Fig. 1. Examples of FORS polarimetric exposures of RS Pup, 
each corresponding to one of the four 60 s exposures labeled 
#12 and #16 in Table [T] Each horizontal stripe corresponds al- 
ternatively to the ordinary and extraordinary beams. The CCD 
gap is visible below the center of the field. The field of view is 
4.3'x4.3^ 



Spectrograph FORSQ (Appenzeller et al. fT998l O'Brien [2009]). 
This instrument is the result of the merging of the original 
FORSl and F0RS2 instruments, and is presently installed at 
the Cassegrain focus of the Unit Telescope 1 (UTl) of the Very 
Large Telescope (VLT). The imaging polarimetry with this in- 
strument is obtained by inserting in the light path a half-wave 
superachromatic retarder plate (hereafter HWP) and a Wollaston 
prism. The HWP introduces a 180° phase shift on the compo- 
nent of the incoming polarized light which is along its reference 
axis, leaving the normal component unchanged (e.g. Clarke & 
Grainger 1971 ). This eff'ectively rotates the incoming polariza- 
tion position angle (PA) by an amount twice as large as the an- 
gle between the incoming PA and the reference axis PA. This 
is equivalent to rotating the whole instrument, and hence the 
Wollaston, by a corresponding amount. The HWP is mounted on 
a rotating support. The Wollaston prism acts as a beam- splitting 
analyser, i.e. it splits the light in two orthogonally polarized 
beams propagating in two directions, separated by an angle of 
22'' on the sky. A strip mask is inserted in the entrance focal 
plane and then re-imaged on the image plane detector, to avoid 
that the two beams overlap. A single exposure therefore gives 
a 50% field coverage only, but with the two polarization direc- 
tions imaged simultaneously on the detector. This optical setup 
including a Wollaston has the advantage of being insensitive to 
changes in the atmospheric turbulence or sky transparency, as the 
two orthogonally polarized images are obtained simultaneously 
(Fig. [T]). We covered the RS Pup nebula with two sets of expo- 
sures (for each epoch), shifted on the sky by 16'' to cover most 
of the nebula. This strategy leaves a series of 6" stripes with- 
out observations, giving a total field coverage of approximately 
80%. The reason for leaving this gap around the central star is 
to avoid its direct imaging on the detector, which would have 
caused a heavy saturation. To cover the rest of the nebula, the 
26 May observation was obtained with the instrument rotated by 
90°. The observations obtained on 26 April are incomplete, with 
only 3 positions of the HWP. This leaves two complete epochs, 
4 April and 26 May, that we will refer to as "first" and "second" 
epochs in the following. 

We obtained four 60 seconds exposures for each position 
of the HWP (0°, 22.5°, 45°, 67.5°), resulting in a total ex- 
posure time of 960s per telescope off'set position, or 1920s 
per complete epoch. The complete list of exposures is given 



# UT date (2009) 



MJD 



pC) AM aQ 



1 

2 
3 
4 
5 
6 
7 
8 

9 

10 

11 

12 
13 
14 
15 
16 
17 
18 
19 



RSPup 

04/04 23: 
05/04 00: 
05/04 00: 
05/04 00: 
05/04 00: 
05/04 00: 
05/04 00: 
05/04 00: 



54:55 
00:49 
06:44 
12:39 
18:54 
24:49 
30:44 
36:39 



26/04 23:49:08 
26/04 23:55:41 
27/04 00:00:33 



26/05 23 
26/05 23 
26/05 23 
26/05 23 
26/05 23 
26/05 23 
26/05 23 
26/05 23 



:16:53 
22:46 
:28:39 
34:32 
40:48 
:46:42 
52:36 
:58:29 



54925.997 
54926.001 
54926.005 
54926.008 
54926.013 
54926.017 
54926.021 
54926.026 

54947.992 
54947.997 
54948.000 

54977.970 
54977.974 
54977.978 
54977.982 
54977.987 
54977.991 
54977.995 
54977.999 



0.1540 
0.1541 
0.1542 
0.1543 
0.1544 
0.1545 
0.1546 
0.1547 

0.6848 
0.6849 
0.6850 

0.4082 
0.4083 
0.4084 
0.4085 
0.4086 
0.4087 
0.4088 
0.4089 



1.33 
1.23 
1.06 
1.05 
1.05 
0.97 
1.15 
1.07 

1.28 
1.45 
1.89 

0.98 
0.73 
0.69 
0.80 
0.90 
0.73 
0.70 
0.73 



1.02 
1.02 
1.02 
1.02 
1.02 
1.02 
1.02 
1.03 

1.06 
1.07 
1.07 

1.24 
1.26 
1.28 
1.30 
1.32 
1.35 
1.38 
1.40 



0.0 

22.5 
45.0 
67.5 
0.0 
22.5 
45.0 
67.5 

0.0 

22.5 
45.0 

0.0 

22.5 
45.0 
67.5 
0.0 
22.5 
45.0 
67.5 





Alphard 










20 


26/05 23:59:12 


54977.999 


0.72 


1.41 


67.5 


21 


27/05 00:13:59 


54978.010 


0.87 


1.26 


0.0 


22 


27/05 00:15:29 


54978.011 


0.83 


1.27 


22.5 


23 


27/05 00:16:59 


54978.012 


0.77 


1.27 


45.0 


24 


27/05 00:18:29 


54978.013 


0.79 


1.28 


67.5 


25 


27/05 00:20:21 


54978.014 


0.75 


1.29 


0.0 


26 


27/05 00:21:51 


54978.015 


0.78 


1.29 


22.5 


27 


27/05 00:23:21 


54978.016 


0.83 


1.30 


45.0 


28 


27/05 00:24:51 


54978.017 


0.73 


1.30 


67.5 



Notes. MJD is the average modified Julian date of the exposures, and a 
is the position angle of the HWP. The seeing p in the visible measured 
by the observatory DIMM is also listed, as well as the observation air- 
mass (AM). 



in Table [T] The phases were computed using the reference 
epoch and period Hsted in Paper I: To = MJD 54090.836 and 
P = 41.4389 days. We used the Vhigh filter #114 of FORS, that 
is a broader version of the standard Bessell V band filter (ef- 
fective wavelength of 557 nm, FWHM of 123.5 nm). The pixel 
scale was 0.252"/binned pixel, and the detector was read in the 
standard mode (2x2 pixels binning). We completed our data set 
with the observation of a PSF calibrator, Alphard (HD 81797, 
HR 3748), in order to measure the scattered light halo in t he in - 
strument and subtract it from the images of RS Pup (Sect. [Z2 ). 
The exposure time was set to 0.6 s to match the ADU counts 
level of RS Pup on the detector. Alphard was imaged at the same 
position as RS Pup on the detector, to obtain the same illumina- 
tion of the array. 



2.2. Data preprocessing 

The individual raw images were pre-processed in a standard 
way (bias subtraction, flat-fielding, bad pixel masking) using the 
IRAIJ^ software package. The flat-field was obtained without the 
Wollaston prism. As discussed by Patat & Romaniello (2006), 
the flat-fielding uncertainties are mitigated by the introduction 



http ://w w w. e so . org/sci/f acilities/paranal/instruments/f or s/ 



^ IRAF is distributed by the NOAO, which are operated by the 
Association of Universities for Research in Astronomy, Inc., under co- 
operative agreement with the National Science Foundation. 
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of some redundancy in the polarimetric measurements (we used 
four HWP positions). We checked that the rotation of the HWP 
does not induce any measurable displacement of the image on 
the detector. A significant diff'erential geometric distorsion how- 
ever exists between the stripes corresponding to the ordinary 
and extraordinary beams. The diff'erential image displacement 
reaches approximately 0.25'' (1 pixel) horizontally, and 1.25'' (5 
pixels) vertically. To compensate for this image distorsion, we 
mapped the ordinary and extraordinary beam images to the same 
fiducial image from our 2008 observations with the NTT/EMMI 
instrument (Paper I). As astrometric references, we selected a 
dozen stars over the image, and we applied a polynomial geo- 
metrical transformation using the IRAF wregister command. 
The resulting RMS dispersion of the star positions relative to the 
transformed world coordinate system (WCS) is approximately 
0.12", or half a camera pixel, which is acceptable. 

The observations of Alphard were pre-processed in the same 
way as RS Pup's. Thanks to the very short exposures, very few 
background stars are present in the images. An exponential 
model of the instrumental/atmospheric scattered light halo was 
derived from the Alphard images using ring median estimates of 
the radial dependence of the halo intensity. This modeling is nec- 
essary to avoid the contamination of the RS Pup images by the 
diff'raction spikes present in the Alphard images at a diff'erent po- 
sition angle. We used the halo model to generate synthetic halo 
images for the ordinary and extraordinary beams. They were 
then scaled to the brightness of RS Pup at the epoch of each ob- 
servation, and subtracted from the added beam images (ordinary 
plus extraordinary stripes) before the computation of the Stokes 
parameters (Sect.|2.3|). 



2.3. Polarimetric quantities 

Following Goldman et al. K009K we obtained the degree and po- 
sition angle of linear polarization by calculating the normalized 
Stokes parameters pq = Q/I and pu = U/I (see also Patat & 
Romaniello 20061 and Zapatero Osorio et al. 120051) . They were 
derived from the processed FORS images using the following 
expressions: 

fo(a) - fE(a) 



F(a) = 



0(a) ^ Ma) 



PQ = Q/I =1/2 [F(0°) - F(45°)] 
p^ = u/I = 1/2 [F(22.5°) - F(67.5°)] 



PL 



^4 



pI^pI 



(1) 

(2) 
(3) 

(4) 



where pi^ is the degree of linear polarization, and /o(o^) and fE(a) 
are respectively the ordinary and extraordinary fluxes obtained 
at the retarder angle a. In this expression, p^ carries a bias due 
to the fact that it is a positive definite quantity. It was therefore 
debiased using Table Al from Sparks & Axon (11999b . A de- 
tailed discussion on this particular question was also presented 
by Simmons & Stewart (1985). The statistical uncertainties of 
the measurements were estimated from the measured standard 
deviation cr(F) of the F(a) values derived from the four expo- 
sures obtained at each position angle of the HWP. The propaga- 
tion of the F(a) errors for the degree of linear polarization p^ 
gives the variance cr^(ph)- 



pt cr^(pQ) + pI cr^(pu) 

^ (/^l) = — ^—-2 

Pq^Pu 



The values of cr^ for each pixel in the pg and pu images were 
derived from the variance of each pixel values in our sequence of 
four F images obtained at the same HWP position angle. As the 
accuracy of this method is limited by the small number of F ex- 
posures, there is a relatively large statistical dispersion on the cr^ 
values. We therefore chose to average the cr^{p^) for each pixel 
over its surrounding 3x3 pixel box to derive more accurate esti- 
mates of the actual variance. We checked that the actual variance 
of the pe pixel values over the 9 pixel box is in agreement with 
this computed average value, giving us confidence in the result- 
ing error bars. The polarization angle 6 and its associated vari- 
ance cr^iO) (for statistically significant polarization degrees) are 
derived using the following expressions (Landi Degl'Innocenti 
etal. 120071) : 



- sign(/7f/) arccos 



(T\e) 



(t\pu) 



\pq + pI/pq\ 




(pI^pI)'' 



(6) 



(7) 



(5) 



The uncertainties we obtain are (t(/7l) ~ 1 to 5% and cr(0) ^ 
0.03 to 0.10 radians (2 to 6°) within 1' of RS Pup (sections of the 
nebula with an average surface brightness). 

As pointed out by Patat & Romaniello (120061) . FORS is af- 
fected by instrumental polarization close to the edge of the detec- 
tor. However, the nebula surrounding RS Pup is contained within 
the central 3x3 arcminutes (including its faint extensions), where 
the instrumental polarization is limited to approximately 0.5%. 
Considering the small amplitude of this bias, we neglect it in our 
analysis. 

2.4. Combination of epochs 1 and 2 

We precisely co-aligned our two observation epochs using 
the background stars as astrometric references. We considered 
epoch 2 as the reference map, and we used epoch 1 only to fill 
the missing parts of the p^ and maps. A moving average scal- 
ing coefficient was applied to the p^, 6 and intensity maps of 
epoch 1 so that they match the local average value of epoch 2. 
The moving average normalization was computed over a square 
area of 60 pixels (15"), and applied multiplicatively to each pixel 
of the epoch 1 image. This procedure ensures a smooth transition 
of the two images and reduces the artefacts. We therefore did not 
average the values obtained at the two epochs for each point on 
the nebula, even where they overlap. There are three reasons for 
this choice: 

- The subtraction of the instrumental PSF wing halo is signif- 
icantly more efficient for epoch 2 as for epoch 1, as the PSF 
calibrator image was obtained immediately after the epoch 
2 on RS Pup, under identical seeing conditions, and with the 
instrument rotator positioned the same way. The subtraction 
of the PSF wings for epoch 1 (obtained 7 weeks earlier) was 
done using this same calibrator observation, therefore less 
efl&ciently. 

- Epoch 2 is of significantly better overall quality than epoch 
1 , due mainly to a lower sky background, and also to slightly 
better seeing conditions (see Table [T]). The mismatch in pho- 
tometric depth and in angular resolution between the two 
epochs would result in artefacts on the combined p^ and 
map. 
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Fig. 2. Left: linear polarization degree /?l measured on the nebula of RS Pup. The dashed lines delimit the sector used to compute 
the profile presented in Fig. [3] Right: position angle 6 of the polarization vectors (in degrees), with contour plots from -90 to +90° 
(steps of 30°). The p\^ and 6 maps result from the combination of FORS epoch 2 (reference) and epoch 1 (to fill the missing vertical 
bands). The field of view is 4' x 4', with North up and East to the left in both cases. 



- The averaging of /?l and 6 values is complicated by their pe- 
culiar statistical properties, and the potential resulting biases 
would be difficult to assess. 

The resulting pi^ and 6 maps are presented in Fig.[2j and the map 
of the Signal to Noise Ratio (SNR) Pl/o'(Pl) is shown in the 
left panel of Fig.|4] The typical SNR on the intermediate surface 
brightness parts of the nebula is between 5 and 10 per pixel, and 
reaches more than 50 per pixel on th e ridge-like structure located 
south of RS Pup (see also Sect. |4.4[ ). A representative profile of 
/7l is presented in Fig. [3] This radial curve was computed by 
averaging the measured pi^ values between position angles 160° 
and 200° (measured from North). The standard deviation of the 
PL values for each radius over this sector is shown with dashed 
lines. 

It is interesting to remark that the position angle of the polar- 
ization vector is always perpendicular to the dust-star direction 
over the whole surface of the nebula. Such a behavior is caused 
by the simple geometry of the system, where the central Cepheid 
illuminates very small dust particles from a single point. With 
no particular microscopic orientation of the particles (that could 
be caused by magnetic fields for example), this creates purely 
Rayleigh scattering in central symmetry. As the position angle 
of the polarization vector does not constrain the nebula geome- 
try, we will not use it for the rest of our analysis. 

The diff'erence between the pi^ maps at the two epochs 
(epoch 2 - epoch 1) is presented in Fig. |4] (right panel). The lin- 
ear intensity scale of this image gives the difference between the 
two epochs for each point in the nebula, expressed in number of 
times the combined standard deviation. The average difference 
between the two images is -\-OAcr. The overestimation of the po- 
larization degree at epoch 1 in the southwest quadrant of the 
nebula (red region in the right panel of Fig.|4]) and its underesti- 






40 50 60 

Radius (arcseconds) 



Fig. 3. Radial profile of the polarization degree pi^ towards the 
south of RS Pup, as a function of the angular radius from the star. 



mation along a southeast-northwest axis (blue regions) is likely 
due to the less efficient asymmetric PSF halo subtraction for this 
epoch (exposures #1-8). The slightly higher polarization degree 
found at epoch 2 over the rest of the nebula is most probably the 
result of the better seeing conditions and the significantly lower 
sky background. 

The combined intensity image from epochs 1 and 2 still con- 
tains small square unobserved patches, at the locations where 
the FORS stripe masks of the two epochs cross each other. To 
retrieve the intensity of the nebula on these patches, we used an 
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10 15 20 25 30 35+0 -3 



Fig. 4. Left: Map of the Signal to Noise Ratio (SNR) of the lin- 
ear polarization degree pi^. The intensity scale is logarithmic. 
Right: Difference in measured polarization degree /7L(epoch2) - 
/7L(epochl), normalized to the combined local standard devia- 
tion of the two images. The field of view is 4' x 4', with North 
up and East to the left. 




2500 3000 



Fig. 5. Total scattered light intensity (in FORS ADUs) from the 
combination of FORS epochs 1 and 2, as well as an EMMI V 
band image to cover the remaining gaps (see text for details). 
The field of view is 4.3' x 4.3'. 



EMMI observation obtained on 27 January 2007 in the V band 
(Fig. [T4j see also Paper I for details). The phase of this image 
is not identical to the phases of the FORS epochs 1 and 2, but 
was obtained at phase 0.877, i.e. approximately in the middle of 
the two FORS phases. As the epochs do not all correspond to 
the same phase of the Cepheid, the intensity distribution of the 
light echoes on the nebula is not strictly identical. This results 
in a slight mismatch of some sections of the nebula in intensity, 
and to discontinuities in the combined intensity image. However, 
these discontinuities remain at a reasonably low level and do not 
aff'ect significantly the aspect of the nebula. The combined inten- 
sity image is presented in Fig.|5] 




Fig. 6. Interpolated map of the inear polarization degree pi^ over 
the RS Pup nebula (from Fig. [2j left panel). The field of view is 
4.3' X 4.3'. 



2.5. Interpolation ofthepL map 

The accurate retrieval of the 3D structure and mass of the nebula 
requires that the pi^ map has a coverage as continuous as possible 
of the nebula. Due to the low surface brightness of the outer 
extensions of the nebula, as well as the presence of gaps due to 
our observing strategy, we therefore need to interpolate the pi^ 
values over the nebula. The missing values were computed using 
a moving average over a 60 pixels square box. The pixels where 
/7l is available within this box were weighted by the square of 
their distances to the interpolated pixel. The field stars visible 
in the image were masked before this operation to avoid biases 
on the degree of linear polarization due to the star haloes. The 
combined pi^ image was then median filtered over a 3 x 3 box to 
reduce the apparent pixel-to-pixel noise, and the resulting map 
is presented in Fig. [6] 

3. Light echo contrast from NTT/EMMI observations 

The observations discussed in this Section were already pre- 
sented in details in Paper I. The processed image data set con- 
sists of a series of six B band images of the nebula from which 
the instrumental/atmospheric scattered light halo of the central 
star was subtracted carefully using a PSF reference. These im- 
ages (Fig. 6 in Paper I) show clearly the propagation of the light 
echoes over the apparent surface of the nebula. 

The fit of the photometric variation curve of RS Pup was ob- 
tained by adjusting the phase (p(a, S), mean value l(a, S) and am- 
plitude AI(a, S) of the Cepheid photometric variation curve to 
the flux measured at each position (a, 6) in the image (where the 
SNR is sufficient). We used a classical ;^^ minimization through a 
Levenberg-Marquardt algorithm. The fit was obtained separately 
on each pixel of 5 x 5 pixel boxes. The standard deviation of the 
measurements over the 5 x 5 pixel box, added quadratically to 
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the statistical fitting error bars, was taken as its associated uncer- 
tainty. While the estimation of each phase off'set is more precise, 
this conservative approach is intended to account for possible 
photometric contaminations by other nearby features. 

As a result, we obtain a map of the normalized contrast 
C(a, 6) of the photometric variation over the nebula, relative to 
the contrast of the photometric light curve of RS Pup in the B 
band Cq: 



C(a, S) = — 



(8) 



with Co = A/rspup//rspup = 1.69. The resulting normalized con- 
trast map is presented in Fig Jt] and an azimuthally averaged ra- 
dial profile is shown in Fig. [spXhis map presents rings close to 
the star. The probable reason of their presence is that the phase 
coverage of our EMMI CCD images (Paper I) is not sufficient to 
have a good 3 -parameter fit (phase/amplitude/average value) at 
all locations over the nebula. This is visible in Fig. 8 of Paper I, 
as most of the measurement points (6) were obtained between 
phases 0.5 and 1.0, and only one point is present between phases 
and 0.5. At some locations of the nebula, we therefore do not 
sample properly the maximum and minimum flux phases, which 
makes the fit more difl&cult. In this case, the amplitude is likely 
underestimated, although the phase is probably much less af- 
fected. We use this map in Sect. |4.1| to estimate the thickness of 
the light scattering layer. The phase map will be used in a forth- 
coming article to estimate the geometrical distance of RS Pup. 
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Fig. 7. Map of the light echoes normalized contrast C(a, S) in 
the nebula of RS Pup. The ring-like structure close to the star is 
probably an artefact caused by the limited phase coverage of our 
EMMI images (see text). The field of view is 3' x 3', with North 
up and East to the left. 
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Fig. 8. Azimuthally averaged radial profile of the light echo con- 
trast C on the nebula of RS Pup, as a function of the angular 
radius from the star. The dashed curves show the boundaries of 
the ±(T domain. 

4. Properties of the light-scattering nebula 

4A. Thickness of the scattering layer 

From the observed amplitude of the light echoes C(a,S) 
(Sect. |3]), it is possible to estimate the thickness of the scatter- 
ing layer. A very thin veil of dust will result in a perfect contrast 
of the echoes compared to the star's light curve. A thick dust 
layer will naturally decrease the echoes contrast, as a contin- 
uum of light-scattering particles with diff'erent phases are then 
present on the line of sight, resulting in a smeared light curve 
(Sugerman [20031 Bond & Sparks 12^91). QuaHtatively, the fact 
that we observe contrasted light echoes indicates that the dust 
is spread over a thin surface, compared to the spatial extent of 
a propagating echo c x P (with c the speed of light and P the 
period of the Cepheid), rather than within a thick l ayer. 

Following a similar method as Havlen (119721) , we define a 
range of softened light curves 7^(0, e) for diff'erent dust layer 
thicknesses e by convolving the Cepheid light curve 7(0) (with 
the phase) in the B band with a square function of variable 
width e. The thickness e of the considered homogeneous dust 
layer model is expressed in units of c x 7* = 7180 AU. We define 
the normalized contrast function C(e) as: 



C(e) = — — 

Is(e) Co 



(9) 



where A7^ is the amplitude of the softened photometric light 
curve, Is(e) its average value, and Co = 1.69 the contrast of the 
light curve of RS Pup in the B band. The shape of the C(e) func- 
tion is shown in Fig. [9] Due to the sharp maximum light peak 
in RS Pup's light curve, this function is close to linear. Inverting 
numerically C(e) gives the relation e(C) between the thickness 
of the scattering layer and the observed contrast. 

We computed the light echo contrast C(a,S) from the 3- 
parameter light curve fit on the EMMI image cube (Sect. [3]). We 
consider the parts of the nebula that are relatively close to the 
star, as we aim at small scattering angles. Large scattering an- 
gles, that would occur e.g. close to the edge of a spherical nebula, 
would result in an apparently higher thickness, due to projection 
eff'ects. In order to derive the true thickness of the dust layer in- 
dependently of projection eff'ects, we thus focus our evaluation 
on the nebula within a radius of ^ 30'' from the star. 
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Fig. 9. Model light echo normalized contrast C(e) (see Eq. [9| 
as a function of the equivalent thickness of the scattering layer, 
expressed in units of ex P = 7180 AU. 
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Fig. 10. Dependence of the equivalent geometrical thickness e of 
the light scattering dust layer as a function of the angular radius 
from RS Pup (expressed in thousands of AU). The dashed curves 
indicate the boundaries of the ±cr domain. 



From the observed contrast C(a,S) and model thickness 
e(C), we derived the equivalent geometrical thickness e of the 
scattering layer. Fig. [T0| shows the resulting thickness as a func- 
tion of the angular radius from the star. We obtain a typical 
thickness of the dust layer of ^ ^ 2 + 2 kAU over the consid- 
ered area (within 3C from the star). At the expected distance 
of RS Pup (1.8 + 0.1 kpc, from Fouque et al. 2007, see also ref- 
erences therein), this corresponds to a relative thickness of only 
^1% compared to the typical radius of the visible nebula (that 
extends on > 2\ Sect. 5^). The weakly significant increase of 
the derived thickness with the angular radius from RS Pup could 
be either due to a geometrical projection eff'ect, or to a physical 
property of the nebula. It is however difficult to conclude firmly 
on this question, due to the limited accuracy of the measurement 
and the fact that this curve was computed through an azimuthal 
integration over the irregular nebula surface. 

We conclude that the geometry of the nebula is essentially 
a single, geometrically and optically thin, light- scattering dust 
layer. Moreover, as discussed in Sect. |4.4[ forward scattering is 
far more efficient than backwards scattering. This means that the 



dust layer we observe in our EMMI images is essentially located 
between us and RS Pup, rather than behind the star. In the fol- 
lowing, we thus adopt the model of a thin veil of dust located be- 
tween us and the Cepheid to reproduce our observations. Thanks 
to the high light echo contrast C(a, S) observed over the nebula, 
we neglect the ring artefacts mentioned in Sect. [3] (that are also 
visible in the azimuth averaged curve in Fig.[TO]) in this analysis, 
as the corresponding underestimation of the contrast does not af- 
fect our conclusion that the scattering dust layer is geometrically 
very thin. 

4.2. Polarization model 

The dependence of the linear polarization degree /7l on the scat- 
tering angle 6 may be assumed to be given by the classical po- 
larization phase function for Rayleigh scattering (White J 97^ : 



1 - cos^ 6 



PL 



(10) 



1 + cos^ e 

where /^max is equal to 1 for theoretical Rayleigh scattering, but is 
lower in real astrophysical situations. The Pl(0) function reaches 
its maximum for 6 = 90°, i.e. dust located in the plane of the sky. 
Considering a maximum polarization value of ;?max = 0-50 + 
0.05, Sparks et al. (12008b showed that the Rayleigh scattering 
model is a relatively good match to the observed polarization 
profile of the light echo of V838 Mon. But these authors also 
deduced from their observations an empirical polarization phase 
function (their Fig. 15) that gives a better match to their obser- 
vations and slightly departs from the Rayleigh function (par- 
ticularly for intermediate polarization degrees). The good re- 
producibility of this phase function for the diff'erent observa- 
tion epochs of V838 Mon (i.e. diff'erent echo distances from the 
star) leads us to retain this empirical model for our analysis of 
RS Pup's polarization profile. The dependence of the degree of 
linear polarization as a function of the scattering angle is pre- 
sented in Fig. [TT] Inverting numerically this function gives us 
access to the scattering angle 6, under the assumption that it is 
lower than 90°. This hypothesis is justified by the fact that back- 
ward scattering is much less efficient than forward scattering, 
as shown by the shape of the photometric Henyey-Greenstein 
phase function O(^) represented in Fig. [TT] (see also Sect. |4.4| ). 
As a consequence, we observe essentially the light scattered by 
the dust located between us and RS Pup (as opposed to the dust 
located behind the star). We may observe particularly dense dust 
features located slightly beyond the plane of the sky (i.e. with 
scattering angles > 90°), but the scattered light flux will be 
low compared to forward scattering features. 

4.3. Three-dimensional dust distribution 

Under the assumption that we see only forward scattering dust 
(0 < 90°) it is possible to invert numerically the polarization 
phase function from Sparks et al. (12008b shown in Fig.[TT]to re- 
trieve the scattering angle 6. The overall geometry of the adopted 
single layer scattering model is presented in Fig. [12] The pro- 
jected distance R of each point of the nebula from RS Pup is 
directly measurable on the image, and we can therefore deduce 
Z immediately, knowing 0, from: 



Z(a,6) 



tan [e(a, S)] 



(11) 



In this expression, Z has the same physical unit as R, and it 
can be expressed either as an angle (expressed, e.g. in arcsec- 
onds), or as a linear distance if we consider the distance to 
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Fig. 11. Degree of linear polarization pi^iO) as a function of the 
scattering angle 6 (in degrees). The empirical pi^{6) function 
from Sparks et al. ( 2008 ) is shown as a solid curve (extrapolated 
to small scattering angles), together with a classical Rayleigh 
scattering function with p^^^ = 0.50 (dashed curve). The pho- 
tometric Henyey-Greenstein phase function O(^) is represented 
for g = 0.538 using a dotted curve. 
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Fig. 12. Geometry of the adopted scattering model. 



RS Pup known a priori. We chose this second approach, assum- 
ing a distance of 1.8 ± 0.1 kpc for RS Pup (n = 0.55 ± 0.03 mas, 
from Fouque et al. 120071) . The pixel scale of our FORS polari- 
metric data (0.252'7pix) thus translates into a linear scale of 
2.22 X 10"^ pc/pix = 458 AU/pix at the distance of RS Pup. The 
resulting "altitude" of the dust layer above the plane of the sky at 
RS Pup's distance is shown in Fig. 13 The relative uncertainty 
on this altitude is estimated to ^ 5 to 15% within V of the central 
star. 



As a remark, the true distribution of the dust likely departs 
from our simple, single-layer model at some locations over the 
nebula. This will be the case for instance if there are several 



distinct layers of dust located on the line of sight. In this case, 
our altitude estimate corresponds to the mean of the altitudes 
of the different dust layers, weighted by the product of the dust 
density by the phase function O(^) of each contributing layer. 

4.4. Dust column density 

To derive the scattered light intensity, we adopt the Henyey- 
Greenstein J1941b phase function (Fig.[TT]): 



O(^) 



47r(l-hg2-2gcos6/)3/2 



(12) 



where 6 is the scattering angle (6 = 0° corresponds to forward 
scattering), and g the scattering asymmetry factor for which we 
adopt a value ofg = 0.538 in the V band (Draine [20033; see also 
Draine 2003b and White 1979 ). Eq.[T2]provides us with a way to 
measure the column density of hydrogen atoms n^ (a, S) of the 
scattering material (in H arcsec"^): 



^ ^^ I(a,6) And" 
f CO O(^) cr 



(13) 



where I(a,S) is the average surface brightness of the neb- 
ula as measured at the Earth at celestial coordinates (a, S) (in 
Wm"^yumarcsec^), d is the linear distance from RSPup to the 
scattering material (in meters), F the apparent average flux of 
the star at the Earth (in Wm"^yum), oj the dust albedo, and 
cr is the scattering cross section per H nucleon. We adopt the 
case A PAH/graphite/silicate grain Milky Way model (Ry = 
Av/E(B-V) = 3.1) devel oped by Weingartner & Draine dSUoTI) . 
and described by Draine (l2Q03al) : oj = 0.674 and cr = 3.275 x 
10"^^m^/H in the V bancj^ Considering the dispersion of the 
difl'erent scattering models, we estimate that a reasonable range 
of values for the product 6L)X(Tis±15% (see e.g. Krugel |2QQ9l) . 
This factor impacts linearly the determination of the density of 
the scattering nebula, and we adopt this range as our systematic 
uncertainty due to the chosen dust model. The only remaining 
variable to measure in Eq. 13 to obtain nu is the distance d be- 
tween the central star and the scattering material. We deduce this 
distance from the geometrical model of the nebula presented in 
Fig. 



13 



For this work, we used the EMMI V band image pre- 
sented in Fig. 14 (see also Paper I). We preferred this approach 
to using the FORS image presented in Fig.|2] as the photometric 
calibration of that direct CCD image is more reliable than the po- 
larimetric observations obtained with FORS. The resulting map 



of the column density nu around RS Pup is presented in Fig. 15 



4.5. Dust geometry model validation 



As discussed in Sect. 4.1 we based our reasoning on the hypoth- 
esis that the light scattering dust layer in the nebula of RS Pup is 
optically thin. A test of this hypothesis can be done by com- 
paring the derived column density values nu to the polariza- 
tion degree /?l for each point in the nebula. For an optically 
thin layer and an overall random dust distribution, these two 
quantities should be independent of each other. Fig. [16] shows 
the two-dimensional distribution of the number of pixels within 
each (/7l, ^h) bin. The small number of pixels with polarization 
values below ^3% (horizontal axis) is caused by the mask- 
ing of the central part of the nebula (where the polarization is 
the lowest) that presents artefacts due to the brightness of the 



http://www.astro.princeton.edu/~draine/dust/scat.html 
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Fig. 13. Altitude of the light scattering dust 
layer relative to the plane of the sky at RS Pup's 
distance (assumed to be 1.8 + 0.1 kpc), ex- 
pressed in parsecs, with contour curves spaced 
by 0.1 pc. The field of view is A'xA' , equivalent 
to 2.1 X 2.1 pc at the distance of RSPup, with 
North up and East to the left. The nebular fea- 
tures located further than 1.8' from RS Pup are 
masked (too low SNR). 
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Fig. 14. EMMI image of the nebula of RS Pup in the V band, ob- 
tained at pulsation phase = 0.^11. The intensity scale follows 
a square root law, and is expressed in W m"^ yum"^ sr"^ The field 
of view is 4' x4\ 



central star. Overall, this diagram appears reasonably well pop- 
ulated, which validates our thin layer hypothesis. We note the 
presence of a population of locations on the nebula with simul- 
taneous high polarization degree (^ 30%) and high column den- 
sity (^ 4 X 10^^ Harcsec"^), in the upper right part of the dia- 
gram. They correspond essentially to the ridge area to the south 
of RS Pup, as well as to other dense areas to the northeast of the 
star. 



4.6. The southern ridge 



The maximum densities are reached over the ridge-like structure 
located to the south of RS Pup, at rin = 10^^ H arcsec"^. An en- 
largement of this region of the nebula is presented in Fig. 17 
This roughly linear nebular feature presents relatively high /7l 
values, and is therefore located close to the plane of the sky at 
RS Pup's distance. It is interesting to remark that this does not 
correspond to the photometrically brightest parts of the nebula, 
that reach only ^h ~ 5.0 x 10^^ H arcsec"^. There is also a slight 
spatial shift between the areas of maximum dust density and the 
lowest altitudes relative to the plane of the sky (Fig. 13). The 
ridge structure therefore appears as a dense ribbon of dust, in- 
clined with respect to the plane of the sky. 
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Fig. 15. Column density of hydrogen atoms riu 
of the scattering material around RSPup, in 
unit of 10^^ H arcsec"^, limited to an angular ra- 
dius of 1.8'. The field stars have been removed 
from the image for clarity. The field of view is 
4' X 4', with North up and East to the left. 
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Fig. 16. Map of the distribution of the measured dust density riu 
(vertical axis, linear scale from to 5 x 10^^ Harcsec"^) as a 
function of the measured polarization degree pi^ (horizontal axis, 
linear scale from to 50%). The color coding is proportional to 
the number of pixels in each (/^l, ^h) bin. 



5. Discussion 

5.1. Total mass of the light-scattering material 

As described in Sect. |4.4| we deduced the surface density of the 
dust from the scattered light intensity and the geometry of the 
scattering layer. This provides us with a means to estimate the 
total mass of the dust layer, by integrating over the nebula. We 
consider in our integration only the parts of the nebula that are 
visible in Fig. [15] We prefer not to extrapolate the density at 
the positions where we do not have a measurement, as the dust 
distribution is inhomogeneous, and this could result in a bias. We 
obtain a total number of H nucleons of ^h = 4.1 x 10^^ ± 15% 
within 1.8' from RS Pup (i.e. ^ 200 kAU, or 1 pc at the distance 
of RSPup). 

For their Milky Way dust model with Ry = 3.1, Weingartner 
& Draine (2001) find a total volume per H atom of Vg = 
2.26 X 10"^^m^H"^ for carbonaceous dust and V^ = 3.94 x 
^Q-33 j^3 pj-i £qj. silicate dust (assuming a total carbon abun- 
dance per H nucleus be = 6x10"^, which is their favored value). 
Considering the small relative abundance of carbonaceous dust 
compared to silicates, we neglect its contribution to the total dust 
mass, and we use the dust grain volume V^ in our computations. 
Moreover, the contribution of the very small carbonaceous dust 
grains is essentially observable at infrared and microwave wave- 
lengths, and is negligible in our visible scattered light observa- 
tions. Following Weingartner & Draine (2001 ), we adopt a den- 
sity of 3500 kg m"^ for the silicate dust, which is an intermediate 
value between crystalline forsterite (3210 kg m"^) and fayalite 
(4390kgm-3). 
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Fig. 17. Enlargement of the ridge nebular feature located 
south of RSPup. Top row: EMMI intensity image (left, 
Wm"^yum"^ sr"\ linear scale) and pi^ map (right, with 2.5% 
spacing contours starting at 0). Bottom row: altitude of the dust 
layer (left, in pc, with 0.05 pc contours) and dust column den- 
sity (right, lO^'^Harcsec"^, with 10 unit contours). The field of 
view is 44 x 44'', corresponding to 0.4 x 0.4 pc at the distance of 
RS Pup. 



Based on our measurement of the total number of H nucle- 
ons, we obtain a total dust volume of Vdust = 1 .62x 10^^ m^ , and a 
total Hght-scattering dust mass of Mdust = (5.7 ± 1.7) x 10^^ kg = 
2.9 ± 0.9 Mq. Taking into account the modeling uncertainties, 
we adopt a systematic uncertainty on this value of ±30%. We 
can now evaluate the total gas and dust mass of the circum- 
stellar envelope of RS Pup by assuming the dust-to-gas ratio 
of Mdust/Mgas = 1% deri ved by Draine & Li (2007) for the 
Milky Way. Draine et al. (12007b compared this model to the 
dust observed in a sample of galaxies, and from their results 
we evaluate the uncertainty on this parameter to ^ 30%. We 
therefore estimate the total mass of the nebula within 1.8' of 
RS Pup to Mgas+dust = 290 M© ± 40%, i.e. approximately be- 
tween 180 and 420 M©. As a remark, it is possible to compute 
directly Mgas+dust = ^H X ^H with mn = 1.672 x 10"^^ kg, that 
gives a value of 345 M©. With a dust to gas ratio of 1%, we 
obtain a dust mass of 3.5 M©, in agreement with our value of 
Mdust = 2.9 ± 0.9 M© derived above. 

The total mass Mgas+dust = 290 M© ± 40% should be consid- 
ered a lower limit for the true mass of the interstellar cloud in 
which RSPup is embedded, for two reasons: 1) we do not in- 
clude the faint extensions of the nebula that are located farther 
than 1.8' from the star (see Sect. |5.2| ), and 2) we are sensitive 
only to the dust that is located between us and RS Pup. The mate- 
rial located behind the plane of the sky at the distance of RS Pup 
remains invisible to us, due to the low efficiency of backwards 
scattering (Fig.[TT]). Under the hypothesis that the dust layer has 
a symmetric counterpart behind RS Pup, the true total mass of 



the cloud could be the double of the above figures, and signif- 
icantly more if we include the faint extensions located beyond 
1.8'. 

An estimate of the dust mass around RS Pup was computed 
by Havlen (11972b . He obtains a dust mass of Mdust = 0-4 M© 
based on a simple multiple shell model, and a second value 
Mdust = 0-03 M© using a more indirect approach (as outlined 
by Lvnds [T968]) . These estimates are both lower than ours , but 
his first value is relatively comparable. As shown in Fig. |A.1[ 
many dust features are difficult to observe at relatively large 
separations on the original photographic plates, and this prob- 
ably caused an underestimation of the total mass by Havlen. 
Moreover, faint parts of the nebula, not counted in Havlen 's es- 
timate, exhibit a particularly high dust density (e.g. the ridge lo- 
cated south of RS Pup). 



5.2. Shape of the nebula 

The light- scattering material surrounding RS Pup appears to be 
spread over an irregular surface (Fig. [T5]). This geometry does 
not present a well defined central symmetry relative to the 
Cepheid. Its density is inhomogeneous, with high dust densities 
over some sections of the nebula (e.g. along the ridge located 
south of RS Pup), and mostly void areas (e.g. in the northwestern 
part). We present in Fig. [18] the result of the direct co-addition 
of our FORS images (#1 to 8 and 12 to 19 in Table [T]), each 
exposed 4 x 90 s, giving a total exposure time of 1.6 h. Faint ex- 
tensions are detectable up to at least 3' from the star, and it is not 
excluded that the nebula extends even further away. The back- 
ground appears slightly non uniform due to the intensity scale 
that is chosen to enhance the faintest extents of the nebula. The 
overall geometry of the nebula appears to have two preferential 
axes, towards the northeast- southwest and northwest- southeast 
of the central star. 
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Fig. 18. Combined FORS V band image of RS Pup. The field of 
view is 6' X 6', with North up and East to the left. 



11 



Kervella et al.: 3D structure and mass of the circumstellar nebula of RS Pup 



As proposed in Paper I, the nebula has likely been shaped 
by the past evolution of RS Pup in a pre-existing, higher den- 
sity interstellar cloud. About ^10 to 20Myr in the past, RS Pup 
went through a phase during which it was a massive, most prob- 
ably fast rotating B type dwarf. During this period, the strong 
equatorial and polar winds of the star combined with its ener- 
getic photons to sweep the interstellar dust away from the star. 
This effect may have compressed the dust into the thin veil that 
we observe today. The current very high intrinsic brightness of 
the Cepheid probably also contributes to the shaping of the neb- 
ula, but at a slower rate. The closest dust to RS Pup that we ob- 
serve in our FORS images is located at linear distances of ^ 1 
to 2 pc from the star, which is already relatively far away. The 
fact that we do not observe a physical displacement of the neb- 
ular features over several decades (Appendix |A]) is consistent 
with this large separation between the star and the dust features. 
Such cavities carved by winds and/or radiation pressure have 
been studied for massive hot stars and supernova remnants (see 
e.g. Lozinskaia ll992K but not so much for regular giant stars. 
It is interesting to remark that the creation of cavities by stellar 
winds and radiation create hydrodynamic instabilities and turbu- 
lence that seem to be present on the surface of RS Pup's nebula 
as wisps and curls. 

Another important phase of RS Pup's evolution was when 
it was a red supergiant (RSG), after its first crossing of the in- 
stability strip and before entering again the instability strip as a 
Cepheid. A plausible scenario could therefore be that mass-loss 
took place during the RSG phase of RS Pup, leading to the for- 
mation of dust relatively close to the star. A good match of the 
position of RSPup in the HR diagram a few 100 000 years in 
the past is provided by Betelgeuse, as the mass of this RSG, 
although uncertain, is comparable to that of RS Pup: Neilson 
& Lester (2011) recently proposed a mass of 1 1.6^3-^ M© for 
Betelgeuse (but Dolan et al. 12008 l obtained ^21 M©), to be com- 
pared with 12.8 ± 0.9 Mo for RSPup (Caputo et al. 2005). The 
subsequent re-heating of the stellar surface after the RSG phase 
would have resulted in the disappearance of the dust close to the 
star, and the sweeping of the interstellar medium. Considering 
the thinness (Sect. |4.ll ) of the observed dust layer, the hypothesis 
of a sweeping of the interstellar dust during the main sequence 
and/or RSG phases of RS Pup appears compatible with the ob- 
servations. 



5.3. Origin of the nebula 

The high mass derived for the nebula is clearly incompatible 
with the scenario that RS Pup created its circumstellar nebula 
through mass loss. The bulk of the material observed in scattered 
light therefore appears to be of interstellar origin. The higher 
density environment in which RS Pup is presently located could 
either be the remnant of the molecular cloud from which the 
star formed, or unrelated interstellar material into which RS Pup 
would be temporarily embedded, due to its proper motion in the 
Galaxy. On the la tter scenario, it is interesting to notice that 
Westerlund (11963b pointed out that the reflection nebula sur- 
rounding RS Pup is a part of a stellar association consisting of 
early-type stars. Although the physical presence of RS Pup in 
this association is possible, the very young age of the associa- 
tion (^4 Myr) compared to the Cepheid may indicate that the 
Cepheid did not form within this association. A more detailed 
discussion on this question can be found in Paper I. 

However, a fraction of the nebular material could still orig- 
inate in past or present stellar mass-loss from RS Pup. As dis- 
cussed previously, mass-loss probably happened during the RSG 



phase that took plac e in th e recent evolution of the star. In addi- 
tion, Kervella et al. ('2009') identified a compact and hot circum- 
stellar envelope around RS Pup, that is likely due to presently on- 
going mass loss. This also supports the findings of Deasy (11988b , 
who identified a very significant mass loss rate for RS Pup (M = 
10"^ Mq yr"^), based on lUE spectra. 

The strong changes in the pulsation period of RSPup 
(Berdnikov et al. 2009 ) might also indicate presently ongoing 
mass loss. The O - C diagram of RS Pup (Fig. 1 in Berdnikov 
et al. 12009b indicates sudden and strong fluctuations super- 
imposed on the secular increase in the pulsation period. The 
secular variations indicate period lengthening at a rate of 137 
s/year, while the most recently superimposed fluctuation in the 
pulsation period corresponds to a sudden change by 0.15% at 
about JD 2 452 000. The recent adaptive optics observations by 
Gallenne et al. (12011b also indicate that the close environment 
of RS Pup contains a significant quantity of hydrogen. The ob- 
served contribution of this envelope in the near-infrared hydro- 
gen lines is significant, and points at a circumstellar location. 

There is another piece of evidence that might refer to the 
origin of some of the material surrounding RS Pup. The radial 
velocity observations available in the literature indicate varia- 
tions in the center-of-mass radial velocity of RS Pup as if this 
Cepheid belonged to a spectroscopic binary system. The com- 
panion could trigger an enhancement of the mass loss from the 
low gravity pulsating atmosphere of the supergiant Cepheid vari- 
able. This companion must have a significantly lower mass than 
RS Pup, as the flux contribution from a more massive star would 
either be easily detectable, or it would have already exploded 
as a supernova (and the nebula of RS Pup does not present the 
characteristics of a supernova remnant). 

6. Conclusion 

From a combination of polarimetric and direct imaging of 
RS Pup's nebula, we derived the geometry of the light- scattering 
dust layer surrounding the star, as well as a map of its column 
density. The observed distribution of the dust appears irregular 
and not centered on the Cepheid. The dust mass of the neb- 
ula is estimated to Mdust = 2.9 ± 0.9 M©, and its total mass to 
Mgas+dust = 290 ± 1 20 Mq. These high values exclude the pos- 
sibility that the bulk of the material of the nebula was ejected 
by the Cepheid itself. The lack of central symmetry in the neb- 
ula also points at an interstellar origin for the nebular material. 
Following Kervella et al. (12009b , we propose that the nebula 
is either a remnant of the interstellar cloud in which RS Pup 
formed, or a particularly dense interstellar cloud in which it is 
temporarily embedded due to its proper motion. There are in- 
dications in the nebula geometry that the Cepheid participated 
in the shaping of the dust cloud, e.g. through its radiative pres- 
sure or wind. The Cepheid may also have contributed a fraction 
of the material present in the nebula, in particular during a past 
RSG phase, but this fraction must be small compared to the in- 
terstellar material, considering the involved dust and gas masses. 
The presence of the large dusty circumstellar nebula around 
RS Pup therefore does not appear to be in itself a property of the 
star, in the sense that it is not the result of stellar mass-loss by the 
Cepheid. Instead, the position of RS Pup in the nebula appears 
essentially coincidental, thus providing a natural explanation to 
the scarcity of similar Cepheid-nebula associations. Although 
bow shocks have been observed by Marengo et al. (2010) in the 
S Cep system, and infrared excess has been identified by Barmby 
et al. (12011b in several Cepheids, the nature of the envelope of 
RS Pup appears fundamentally diff'erent, both in extension and 
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total mass. Another consequence is that the dusty nebula is not a 
good candidate to resolve a possible evolutionary vs. pulsational 
mass discrepancy for RSPup (see e.g. Keller 2008, Caputo et 
al. 12005, and Neilson & Lester 2008, 2009). However, the com- 
pact and hot envelope observed by Kervella et al. (,2009J still 
appears as a plausible contributor. 

Our next step, that will be the subject of a forthcoming arti- 
cle, is to determine the distance of RS Pup using its light echoes, 
as the geometry of the light-scattering dust is now well deter- 
mined. 
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Appendix A: Search for temporal evolution of 
RS Pup's nebular features 

A.1. Westerlund's photographic plates 1960-1963 

Soon after his discovery of the nebula surrounding RSPup, 
Westerlund obtained a series of photographic plates with the 74- 
inch Mount Stromlo telescope (located near Canberra, Au stralia ) 
to search for the propagation of light echoes (Westerlund 11961b . 
These plates were subsequently used by Havlen for his pioneer 
study of the nebula (Havlen 11972b , and a detailed description 
of the plates can be found in this publication. For the present 
work, we digitized the original photographic plates, that were 
preserved by one of us (RJH) for almost 50 years, using a flatbed 
scanner at a resolution of 2400 dpi. This optical resolution is 
sufficient to resolve the film grain. We matched the astromet- 
ric WCS of the scanned plates with the EMMI images obtained 
in 2007 using a dozen of background stars spread over the im- 
age field. The proper motion of a few bright stars of the field is 
noticeable, however, the matching of the images is satisfactory. 
We carefully blinked the Westerlund plates with the EMMI 
observations presented in Paper I. The EMMI images were in- 
terpolated at the phase of the Westerlund observations using the 
procedure described in Sect. [3] Among the 18 epochs cover- 
ing the 1960-1963 period, we focused our comparison on three 
epochs that show particularly well the circumstellar nebula: 30 
March 1962 (visual), 5 May 1962 (visual) and 26 February 1963 
(blue and visual). These plates were obtained under very good 
seeing conditions, and are the best of the series. We did not de- 
tect any change in the morphology of the nebula. This negative 
result was expected considering the fact that a T' angular dis- 
placement corresponds to a linear displacement of ^ 2000 AU 
at the distance of RSPup. Over 50 years, such a V displace- 
ment would therefore require a very fast linear velocity of at 
leasts 200 km/s. 

A.2. CFHT observations 1992-1993 

The circumstellar nebula of RS Pup was observed on the nights 
of 4-6 January 1992, 9 April 1992, and 28 March 1993, using 
the HRC instrument installed at the prime focus of the Canada- 
France-Hawaii Telescope, for Michael J. Pierce & Douglas L. 
Welch. Exposures of 30 to 600 s were obtained in the B, V and 
R bands. We present in Fig. A.2 the best image in the V band 
(06-01-1992), compared to our second epoch observation with 
FORS (27-05-2009). The two observations do not correspond to 
the same phase, but we could not detect any change in the mor- 
phology of the nebula (apart from the photometric variations due 
to the echoes). However, we observe the apparent displacement 
of several stars in the field over the 17 years separating the two 
observations (an example is presented in Fig. ??). The bright 
pair of point-like sources visible above the center of the FORS 
image (close to the central vertical white stripe) is an artefact of 
the instrument. 
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Fig. A.l. Comparison of Westerlund's photographic plate #31 17 with the interpolated and PSF- subtracted EMMI observation at the 
same phase (0 = 0.84). The field of view is 1.4' x 1.4', with North up and East to the left. 




Fig. A.2. CFHT/HRC image of the circumstellar nebula of RS Pup (left, 1992) and the same field observed with FORS (right, 2009). 
The field of view is T x 2', with North up and East to the left. 
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